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Abstract. We propose an empirical sulphur abundance 
determination method based on the strong sulphur emis- 
sion lines, [Sn] AA6716, 6731 and [S in] AA9069, 9531 in 
H II regions. From a compilation of literature data we have 
made a calibration of sulphur abundance versus ( [S II] + 
[Sm])/H/3, similar to what has been widely used for the 
' more easily observable iso-electronic element oxygen. This 
enables abundance determinations in extragalactic H II re- 
gions without measurements of weak temperature sensi- 
tive lines for use in model calculations. 

As a first application of the empirical calibration 15 
spectra covering the wavelength range 3650-10000 A have 
been obtained of H n regions at various galactocentric dis- 
tances in the spiral galaxy NGC 300, by employing spe- 
cially designed multiapcrturc slitmasks which allow simul- 
taneous observations of several Hii regions. Sulphur and 
oxygen abundances are determined and an oxygen gradi- 
ent in agreement with previously published work is found 
as well as a somewhat steeper sulphur abundance gradient. 
This results in a slight decrease of log S/O with increasing 
radius. 

The quite low values of S/O seem to confirm the trend 
demonstrated by Diaz et al. (1991) of a decrease in log 
S/O at high metallicities, but still more data are needed 
before any definite conclusions can be drawn. We find as 
expected that the radiation softness parameter rj decreases 
with increasing radius and the slope is comparable to what 
is found for spiral galaxies of similar metallicity. 

Key words: H n regions - galaxies: abundances - galax- 
ies: NGC 300 - galaxies: spiral 
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1. Introduction 

Elemental abundance studies in spiral galaxies have 
largely been concerned with the abundance of oxygen and 
the presence of radial O-gradients with negative slopes 
have been demonstrated for more than a decade. Initiated 
by the analyses and line ratio calibrations of Edmunds & 
Pagel (1984) and McCall et al. (1985) large observational 
samples have recently been compiled by several groups 
(Vila-Costas k Edmunds 1992; Oey & Kennicutt 1993; 
Zaritsky et al. 1994, hereafter ZKH), of which the latter 
collates much of the existing published and their own spec- 
troscopic data on a total of 577 Hn regions in 39 spirals. 

The present data available through these surveys are 
of such proportions, that they can address the correlation 
of oxygen abundance with local and macroscopic param- 
eters of spiral galaxies, which can be used to constrain 
galaxy evolution models. ZKH define a characteristic oxy- 
gen abundance and find that it correlates with galaxy hub- 
ble type and luminosity. No clear correlation of oxygen 
abundance gradients, when expressed as dex/isophotal ra- 
dius, with Hubble type is seen, although there are indica- 
tions that intermediate-type spiral galaxies have the steep- 
est gradients and that the slope is flatter in barred than 
unbarred spirals. 

Much less attention has been paid to the sulphur con- 
tent of spiral galaxies, probably because the strong lines 
of the ionization stage S 2+ is situated in the near-IR re- 
gion at A = 9069 and 9531 A, where detector sensibility 
is low and many atmospheric emission lines are present, 
although it can give important clues on elemental abun- 
dance and the softness parameter as defined in Section 
5.1, which was shown by ZKH to increase with decreasing 
radius, etc. 

Garnett (1989) investigates sulphur abundance in 13 
extragalactic H n regions and summing up on his own and 
previous work concludes that there is only very weak if any 
evidence for gradients in S/O in spiral galaxies. Neither 
is there a systematic variation of S/O with oxygen abun- 
dance even though it might be expected from chemical 
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evolution models. In this work we will put more emphasis 
on the sulphur abundance in spiral galaxies by developing 
an empirical calibration of sulphur line ratios and abun- 
dance. 

It has become evident from the studies cited above, 
that conclusions on abundance variation over a galaxy 
disk drawn on data from less than five H II regions can 
be very misleading and even turn out with the wrong sign 
of the gradient. Therefore more than 10 observed regions 
per galaxy with good radial coverage is desirable. Ryder 
(1995), who added oxygen data on 6 new galaxies with 
findings in support of ZKH, has discussed how Hii re- 
gions can be observed in a time efficient manner with nar- 
row band filters and multifiber or multislit spectroscopy 
and pointed out some of the drawbacks of the first two. 
When broad wavelength coverage of H n regions in nearby 
galaxies is the objective, multislit spectroscopy might be 
prefered, since most regions extend over many arcseconds 
and may be closely situated. Hence, they are not effi- 
ciently observed using multifiber spectrographs designed 
for point-like objects. 

NGC 300 is a large, nearby spiral galaxy with many 
prominent Hii regions. At the modest distance of the 
Sculptor group of galaxies, about 1.2 Mpc, the Scd galaxy 
extends over 20' and is a perfect object for multislit spec- 
troscopy covering wide fields. By this technique using spe- 
cially designed multiaperture plates the number of H II re- 
gions that can be observed simultaneously is increased as 
compared to conventional long-slit spectroscopy. 

2. Line ratios as abundance indicators 

Abundances from H II region emission lines are most ac- 
curately determined when temperature sensitive lines are 
measured, since the derived temperature enables model- 
ing of the region. The required lines such as [O in] A4363 
are, however, quite weak and often not detectable at all. 
In their absence one has to depend on empirical relation- 
ships between the desired abundance and a combination 
of strong emission lines. Pagel et al. (1979) suggest R23 = 
([On] AA3726, 3729 + [Om] AA4959, 5007)/H/3 as an im- 
proved oxygen abundance indicator as compared to earlier 
used line ratios. 

A calibration of log(R23) against oxygen abundance by 
Edmunds & Pagel (1984) has since been widely used for 
abundance determinations in regions where the tempera- 
ture sensitive lines are unreachable. The fit is ambiguous, 
however, consisting of two branches, and care must be 
taken to make sure whether the high or low mctallicity 
branch is appropriate in each case. The uncertainty of the 
abundance is ±0.2 dex. 

2.1. A sulphur abundance indicator 

The success of this empirical method of oxygen abun- 
dance determination has inspired a similar method for sul- 



phur abundance determination. Since these elements are 
iso-electronic, an expression corresponding to R23 can be 
made for sulphur, simply by replacing each oxygen line by 
the corresponding sulphur line. We define: 

_ [Sn] AA6716, 6731 + [Sm] AA9069, 9531 , , 
S23 = p • (1) 

As two of the sulphur lines arc in the far-red part of 
the spectrum, these have only in the recent years become 
available, as detector sensibility in this spectral range has 
been improved. It might be more relevant to normalize S23 
to Ha, the nearest hydrogen Balmcr line, than to H/3 lying 
further away in the blue range where it is more sensitive 
to extinction correction uncertainties. Our spectral reso- 
lution, however, is so poor that we can hardly deblend Ha 
from the nearby [Nil] A6548 and [Neil] A6583 lines and 
thus we find it more reasonable to use H/3. A calibration of 
S23 versus sulphur abundance should now enable sulphur 
abundance determinations in H 11 regions where the weak 
temperature sensitive line [S ill] A6312 is undetectable. 

A search through literature for the data available for 
such a calibration has been made. The resulting plot is 
shown in Fig. 1. All papers have measured values of emis- 
sion line intensities of the constituent lines of S23 and give 
from ionization model calculations the sulphur abundance. 
The data points from Diaz et al. (1991) are from giant 
Hn regions in M51, probably the spiral with the high- 
est abundance known. Their position in the plot suggests 
that they belong to a high metallicity branch, whereas all 
other data points seem to fit nicely on a low metallicity 
branch. It seems likely that only extremely high metal- 
licity H 11 regions will be on the high metallicity branch. 
The majority of H 11 regions will fit on the lower branch 
and thus only for this branch a fit is made. A first order 
weighted fit gives: 

log(S/H) + 12 = 6.485 + 1.218 log(S 23 ). (2) 

In addition a grid of model H 11 regions from Stasinska 
(1990) has been plotted in the diagram. All models have 
one ionizing star and both models with uniform hydrogen 
number density 10 and 1000 cm~ 3 have been included. 
The horizontal dispersion is mainly due to the effective 
temperature of the ionizing star varying from 32500 to 
55000 K. The theoretical models further consolidate the 
lower branch and suggest the existence of a higher branch. 
The latter, however, has as yet a quite diffuse appearance. 

All papers used for this plot except Pastoriza et al. 
(1993) have taken into account the possible contribution 
of S 3+ to the sulphur abundance. All but one of their 
regions have so large + /0 ratios (> 0.5) that S 3+ is not 
likely to be present in significant amounts (Garnett 1989). 
The rightmost of the datapoints in Fig. 1 has considerably 
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Fig. 1. Sulphur abundance as a function ofS23- x : Diaz 
et al. (1991). * : Garnett (1989). A : Garnett & Kennicutt 
(1994). O : Osterbrock et al. (1992). + : Pastoriza et al. (1993). 
□ : Skillman et al. (1994). • : Theoretical models from Stasin- 
ska (1990), not included in the fit. The dashed line is the fit 
given in Eq. (2). 

higher ionization (as it is a region in the nucleus) , and S + 
+ S 2+ might not contribute more than 80 % of the total 
sulphur abundance. 

3. Observations 

The observations were carried out at the Danish 1.5-m 
telescope at ESO, La Silla, Chile, where a facility to ex- 
ploit multislit spectroscopy with the Danish Faint Ob- 
ject Spectrograph and Camera (DFOSC) is under devel- 
opment. The CCD available at the time of observation 
(Thompson lk x lk) had a field size of 8.'6 x 8'6, large 
enough to include the main parts of the galaxy. The aper- 
ture plates are manufactured on site from a hard plas- 
tic material with a milling machine controlled by an in- 
struction file generated on the basis of well-sampled Ha- 
exposures obtained with the same instrument. 

NGC 300 was observed during 2 nights in July 1995 
(see Table 1). Two different aperture plates were used, 
both having 8 non-overlapping slitlcts with widths of 
~ 2'.'5. The number and lengths of the slitlets were gov- 
erned by the wish to maximize the number of Hn re- 
gions observed in one exposure and the need to have suf- 
ficiently sampling of the galaxy stellar continuum and the 
local sky background, which is especially important in the 
near-IR range where emission from atmospheric OH-lines 
rises strongly. Since some regions have large angular sizes 
and some slitlets included two regions these criteria were 
meet using lengths of 30-60". One slitlet was positioned 
to contain pure sky background for comparison, but the 
use of such a global sky spectrum gave generally a poorer 
background subtraction, presumably due the low spectral 



resolution (see below) and the displacement of the spectra 
in the dispersion direction. 

In total spectra of 14 H n regions were obtained with 
the two aperture plates; region 137A appears in the ob- 
servations from both nights due to an initial misidentifi- 
cation. With exposure times of 3000-7000 s the advantage 
of simultaneous observation of 5-10 Hil regions compared 
to long-slit spectroscopy is obvious. Unfortunately the 
achieved spectral resolution was rather low, ~ 25 A and 
<~ 60 A in the blue and red ranges respectively, due to 
a malfunction of the best suited grisms. This means that 
[S II] AA6716, 6731 must be deblended from the Hci line 
at 6678 A and the flux of weak Paschcn lines must be 
subtracted from the [S in] lines. 

4. Data reductions 

The multiaperture CCD frames were reduced by employ- 
ing standard IRAF packages for flat-fielding, extraction of 
one-dimensional spectra, calibration etc. Since the slitlets 
cannot be milled with exactly the same widths special care 
must be taken in the reduction process. Pixel-to-pixel vari- 
ations are removed by normalizing the spectrum from each 
slitlet using well-exposed dome flat-fields made through 
the same instrumental setup. Long range variations over 
the CCD are accounted for by normalizing the collapsed 
one-dimensionel spectra from each slitlet of a twilight sky 
exposure with the signal of the slitlet through which the 
flux standard star exposures are taken. This last proce- 
dure is similar to the illumination correction applied in 
long-slit spectroscopy. 

The wavelength calibration was performed on the ex- 
tracted one-dimensional spectra using He lamp exposures, 
both for the blue and red wavelength range, obtained 
through the slitlet mask right after the science frames, 
while the telescope was pointing in the same direction 
to minimize any effect of flexure. Multiple exposures of 
the spectrophotometric standard star LTT 377 (Stone & 
Baldwin 1983, Baldwin & Stone 1984) from each night 
were used to convert to absolute flux. The effect of at- 
mospheric refraction (Fillipcnko 1982) was negligible with 
the employed slit width, since both NGC 300 and LTT 377 
were observed at low airmasses ranging between 1.03-1.08. 
The blue and red spectra overlap in the region 5300-6200 
A, which can be used to test the reliability of the flux cal- 
ibration. Generally the deviation between the two ranges 
is less than 5 % (see Figure 2) and the mean difference for 
the apertures with high S/N is ^ 2% for both nights, so 
no correction factor are applied when the two parts of the 
spectrum are merged. 

In the near-IR spectral range careful measures were 
taken to remove the effect of the atmospheric H 2 ab- 
sorption band at AA8900-9800 using the procedure pro- 
posed by Osterbrock et al. (1990). After fitting a smooth 
global continuum the line fluxes were measured by fitting a 
Gaussian profile to the line, or two Gaussians if deblend- 
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Table 1. Log of observations 



Date 


Hn regionID a 


Wavelengths (A) 


exp. time (s) 


July 30 1995 


45,61,77,79,109,118A,127,137A 


3650-6200 


2700 


July 30 1995 


do. 


5300-10000 


4500 


July 31 1995 


53AB,76A,88+90,100,H9A,137A,137C 


3650-6200 


4800 


July 31 1995 


do. 


5300-10000 


7200 



a Identification of the individual Hn regions follows the numbering by Deharveng et al. (1988) 



ing from nearby lines is relevant, and integrating under 
these profiles. 

Finally the measured line fluxes are corrected for in- 
terstellar extinction by determination of the visual ex- 
tinction, Ay, using the line flux of H<5, H7, H/3 and Ha, 
and where possible the Paschen line P9, with the method 
described by Petersen & Gammelgaard (1996) assuming 
T c = 10 4 K and n c = 10 2 cm~ 3 for the giant H 11 re- 
gions. For those H II regions where the extinction deter- 
mination does not seem to be reliable, values of Ay have 
been adopted from D'Odorico et al. (1983) and Webster 
& Smith (1983) if available, or as a simple average of the 
other regions. 

The data reduction can be tested by comparing the 
ratio of the two near-IR [S in] lines with the theoretical 
ratio of 2.48 (Mendoza & Zeippen 1982). We find a mean 
value for 7(A9531)/7(A9069) of 2.78 ± 0.15, i.e. 12 % above 
the theoretical ratio. In fact this deviation is quite small 
giving us faith that our correction for the atmospheric 
H2O absorption band and blending by OH emission lines 
has been successful in spite of the low spectral resolution 
we were forced to use. This is probably due to the fact 
that both [Sm]A 9069 and [Sill]A 9531 are fortuitously 
positioned at wavelengths where atmospheric effects are 
relatively low as evident from the plots in Osterbrock et al. 
(1990) and Osterbrock & Martel (1992). We are therefore 
confident that the sulphur line fluxes can be assessed with 
sufficient accuracy for the present study. 

5. Results 

An example of the final reduced spectra is shown in Fig. 
2. In the region 5300-6200 A, covered by both the blue 
and the red exposure, the success of the flux calibration 
is demonstrated by the closely coinciding spectral parts. 
In Table 2 the extinction and water absorption corrected 
line intensities relative to H/3 are listed along with the 
absolute H/3-flux and the extinction in V. The underly- 
ing Paschen hydrogen emission has been subtracted from 
the [S ill] line fluxes by means of the theoretical line ratios 
from Hummer & Storey (1987), assuming physical proper- 
ties as mentioned above for the extinction determination. 

An exact uncertainty calculation for the line intensi- 
ties is very hard to carry out due to the many sources 
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Fig. 2. Combined spectrum of the region 53A, B. Notice how 
smoothly the blue and red part of the spectrum merge. In the 
overlapping region the difference in resolution is obvious. The 
most important lines are labeled. 



of uncertainty in the reduction steps and the deblending 
procedure. Consequently we have assigned uncertainties 
to the lines following these qualitative guidelines inspired 
by Ryder (1995): when S/N is high and the lines easily 
identified we assess a relative uncertainty of 15 % and 
when the lines could only just be distinguished from the 
noise we adopt a relative uncertainty of 30 %. In the case 
of Ha the uncertainty assigned is 30 % due to problematic 
deblending from nearby lines. 

The agreement between the observations from the two 
nights might be tested by comparing the two measure- 
ments of region 137A. The spectrum from July 31st (listed 
first in Table 2) has somewhat higher resolution and the 
value of A v derived from these data has been used for 
both spectra. The relative line intensities for [O n] and 
[S n] are lower the second night while the lines from the 
double-ionized stages are considerably higher (20 % for 
[O in] and 35 % for [S in] ) . This is more than the expected 
uncertainty and suggests that the two slitlet positions do 
not coincide as we see a higher ionization the second night 
and furthermore a much lower absolute H/3-intensity. 
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Region: 


45 


53 A, B 


61 


76A 


77 


79 


88, 90 


100 


fO III AA3726 372Q 


3.441 


3.507 


2.553 


3.103 


2.365 


3.018 


4.366 


2.034 a 


R8 


0.293 a 


0.294 










0.419 




ii j 


0.477 


0.494 




0.055 a 










H/3 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 a 


ro ml A4959 


0.105 a 


0.549 


0.481 


0.123 a 


0.464 


0.380 


0.254 


0.221 a 


[Oin] A5007 


0.368 a 


1.580 


1.503 


0.297 a 


1.421 


1.105 


0.696 


0.789 a 


Hq 


1.81 a 


2.25 a 




2.28 a 






2.18 a 




[Sii] AA6716, 6731 


0.623 


0.344 


0.292 


0.478 


0.376 


0.458 


0.578 


0.821 a 


[S m] A9069 


0.168 a 


0.217 


0.214 


0.285 


0.164 


0.175 


0.169 


0.557 a 


P9 


0.042 a 


0.038 




0.038 a 










[S m] A9531 


0.331 


0.484 


0.516 


0.860 


0.348 


0.370 


0.527 


1.807 a 


-log /(H/3) 


13.07 


13.96 


13.62 


14.51 


13.19 


13.28 


13.96 


13.43 


A v 


1.04 


0.91 


0.87 d 


0.89 


0.86 b 


1.01 b 


1.23 


0.79 c 



Region: 


109 


118A 


119A 


127 


137A 


137A 


137C 


[0 n] AA3726, 3729 


3.293 


2.205 


3.151 


1.955 


1.750 


2.094 


2.523 


ns 




0.226 


0.254 a 




0.321 a 






H 7 




0.497 


0.435 




0.429 




0.497 a 


H/3 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


[0 in] A4959 


0.163 a 


0.597 


0.348 


0.747 


0.862 


0.724 


0.625 


[0 in] A5007 


0.588 a 


1.771 


1.041 


2.293 


2.583 


2.125 


2.030 


Hq 






2.80 a 




2.18 a 




2.31 a 


[Sii] AA6716, 6731 


0.582 


0.291 


0.467 


0.302 


0.162 


0.172 


0.224 


[S in] A9069 


o.ior 


0.175 


0.223 


0.220 a 


0.270 


0.198 


0.114 a 


P9 




0.028 a 


0.024 a 




0.027 a 




0.032 a 


[S in] A9531 


0.303 


0.586 


0.845 


0.462 


0.775 


0.578 


0.406 


-log I (H/3) 


13.42 


13.17 


13.74 


13.70 


13.14 


13.79 


13.18 


A v 


1.22 b 


0.79 


0.42 


0.87 d 


0.67 


0.67 


1.70 



a Uncertainty estimated to 30 %. Elsewhere 15 %. 
b From D'Odorico et al. (1983). 
c From Webster and Smith (1983). 
d Average of the other regions. 



From the line intensities the parameter S23 is calcu- 
lated according to Eq. (1) and plotted versus radial dis- 
tance from the center of the galaxy in Fig. 3. The depro- 
jected fractional isophotal radii p/p , where the isopho- 
tal radius is the distance where the surface brightness 
equals 25 mag/arcsec 2 , are taken from Deharveng et al. 
(1988). For NGC 300 p a = 9.'75. This normalization ra- 
dius is less sensitive to selection effects, when comparing 
physical properties of spiral galaxies. The uncertainty cal- 
culation for S23 takes into account the uncertainty of the 
extinction and the estimated uncertainty of the line fluxes. 
From Fig. 3 a gradient in S23 is evident. A linear weighted 
fit gives log S 23 = 0.40 ± 0.09 - (1.03 ± 0.29)p/p o . By 
means of the relation in Eq. (2) this is converted to a sul- 
phur abundance gradient: 



12 + log (J^j = 6.97 ± 0.14 - (1.25 ± 0.38) p/ Po (3) 

For the rather high S23-values we find, one might won- 
der whether the regions could belong to the (as yet in- 
definite) high metallicity branch in Fig. 1 in stead of the 
lower branch where our Eq. (2) is valid. However, since we 
find a decrease in S23 with radius, only the lower branch 
is consistent with a decreasing sulphur abundance gradi- 
ent which ensures us that our assumption is correct. A 
fit to the high branch would be inversely proportional to 
the abundance and would result in a sulphur abundance 
gradient increasing with radius in contradiction with all 
expectations. 

From the observed oxygen line fluxes the parameter 
R23 is calculated. The calibration by ZKH of R23 versus 
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Table 3. Abundances 



Region: 


45 


53A, B 


61 


76A 


77 


79 


88, 90 


100 


P/Po 


0.27 


0.39 


0.32 


0.09 


0.35 


0.35 


0.37 


0.07 


1Z -f- log \ J / 1 1 


ft Q1 

o.y i 


ft 71 


s ft/i 


ft OR 

o.yo 


ft ft7 


ft ft/i 


ft 


Q 09 


1 9 -I- lno- S /H c 
±Z log o / LI 


KJ.OO 


U.0 1 


u.ou 


fi 74 


fi 49 




U.Ul 


7 i n 

1 .1U 


— log 0/ \J 


9 *37 
Z.o I 


9 90 


9 ^4 
Z .O^l 


9 99 
Z . ZZ 


9 

Z .^10 


9 ^ 
z.ou 


o 19 
Z.IO 


1 Q9 
i.yz 


Region: 


109 


118A 


119A 


127 


137A 


137A 


137C 




P/Po 


0.17 


0.29 


0.34 


0.30 


0.40 


0.40 


0.42 




12 + log 0/H b 


8.90 


8.84 


8.84 


8.79 


8.76 


8.79 


8.76 




12 + log S/H c 


6.48 


6.51 


6.71 


6.48 


6.58 


6.46 


6.33 




-log S/O c 


2.42 


2.32 


2.13 


2.31 


2.18 


2.34 


2.44 





a From Deharveng et al. (1988). 

b Uncertainty estimated to ± 0.2 dex. 

c Typical uncertainty: ± 0.3 dex. 



P/Po 



P/Pa 



Fig. 3. The variation of the parameter S23 = ([S n] + [S IIl])/H/3 Fig. 4. The oxygen abundance gradient in NGC 300. The ab- 

across NGC 300. The abscissa gives the radial distance p in scissa gives the radial distance p in units of the isophotal radius 

units of the isophotal radius p a . The dashed line is the best p . The dashed line is the best linear fit given in Table 4. 
linear fit given in Eq. (3) . 



oxygen abundance is employed and the resulting plot of 
oxygen abundance versus radius shows a clear trend of de- 
creasing abundance with growing radius (Fig. 4). In Table 
3 the calculated abundances for each region are given as 
well as the isophotal radii. The calculated abundance gra- 
dients are given in Table 4 together with results for oxygen 
from previous work. There is a fine agreement between the 
results. 

Comparing the gradients in Figs. 3 and 4 it is seen that 
the oxygen abundance gradient is more well-determined 
than the one for sulphur. This is most likely due to the 
much lower resolution of the red spectrum where all the 
sulphur lines are located as compared to the blue spectrum 
where all the oxygen lines lie. 



As the sulphur abundance has a steeper slope than the 
oxygen abundance we find a slight decrease in S/O as the 
radius increases as can be seen in Fig. 5. The value of the 
gradient is given in Table 4. 

From Fig. 3 it is apparent that the value of the sulphur 
abundance gradient is quite dependent on the innermost 
datapoint. If this region (# 100) is omitted the gradient 
is drastically reduced to —0.43 ± 0.37 dex/p , indicating 
that the uncertainty of the gradient might very well be 
greater than that quoted in Table 4. We do not, however, 
especially mistrust our observations of this region. Even 
though it has the weakest emission of our sample the S /N 
of the spectrum is acceptable. To firmly establish the value 
of this gradient observations of more regions at large radial 



Christensen et al.: A Sulphur abundance study of NGC 300 by an empirical calibration method 
Table 4. Abundance gradients in NGC 300 



7 



Source 



Abundance gradient 



This work 
This work 
This work 

Deharveng et al. (1988) 
Zaritsky et al. (1994) a 



12 + log (S/H) = (6.97 ± 0.14) 
12 + log (O/H) = (9.03 ± 0.04) 
log S/O = (-2.12 ± 0.11) 
12 + log (O/H) = (8.95 ± 0.04) 
12 + log (O/H) = (8.97 ± 0.04) 



- (1.25 ± 0.38) p/p 

- (0.63 ± 0.13) p/po 

- (0.50 ± 0.35) p/po 

- (0.63 ± 0.10) p/p 

- (0.61 ± 0.05) p/po 



Based on the combined data from Webster and Smith (1983), Pagel et al. (1979) and Deharveng et al. (1988). 



0.20 0.30 
p/Po 



Fig. 5. log S/O plotted vs. fractional isophotal radius p/p - 
The dashed line represents the best linear fit, see Table 4. 



distances as well as at small distances are needed (prefer- 
ably at a higher resolution). 

The same sensitivity to this innermost datapoint is 
displayed by the gradient in Fig. 5. When the S/O gradient 
is calculated with region 100 omitted the result is +0.12 
± 0.36 dex/,o , i.e. a radial gradient in S/O is no longer 
present. 

If a relation between log S/O and oxygen abundance 
could be established a powerful tool to test galactic chem- 
ical evolution models would be at hand. Sulphur is pro- 
duced in stars with M > 10 Mq and the ratio of sulphur 
to oxygen produced is greatest in stars with masses in 
the range 12-20 solar masses. Variations in log S/O vs. 
metallicity could thus reflect variations in the IMF. Gar- 
nett (1989). Investigations of this relationship has been 
made by Garnett (1989) and Diaz et al. (1991). The for- 
mer find no evidence for a relation while the latter investi- 
gates more metal rich regions and find a decrease at high 
metallicities. Our data confirm the trend demonstrated by 
Diaz et al. (1991) of a decrease in log S/O at high oxygen 
abundance. We do, however, find quite low S/O- values in 
comparison with their collection of data, but keeping the 
rather large uncertainties (including both uncertainties of 



the fluxes and the fit in Eq. (2)) in mind, the discrepancy 
might not be significant. 

5.1. The softness and ionization parameters 

Some of the physical parameters determining the Hll 
region spectra arc the ionization parameter u and the 
ionizing spectrum characterized by the mean effective 
temperature of the ionizing stars. Vflchcz & Pagel 
(1988, hereafter VP) introduced the parameter n = 
(0 + /0 2+ )/(S 2+ /S + ), which was shown to have a one-to- 
one relation with T*, as an indicator of the radiation soft- 
ness. It has the important advantages of being insentive to 
density and virtually independent of u and T . The actual 
relationship with T* depend on stellar model atmosphere 
calculation. 

As discussed by Garnett (1989) it is still premature 
to make any quantitative assessment of from n and 
more detailed non-LTE models of hot stars are required to 
derive absolute values of the effective temperature of the 
ionizing sources. Nevertheless, he argues that r\ can still 
be used as a sequencing parameter to examine the relative 
behaviour of T*, and we shall use it to investigate the gross 
proporties of T± in NGC 300. When the near-IR [S ill] lines 
are available r\ is easily found from the observable line ratio 
quotient (VP): 



77/ = 



[0 11] AA3726, 3729 / [S 11] AA6717, 6731 



[Oin] AA4959,5007 / 
through the relation 
log 77 = 77/ + 0.14/t + 0.16, 



[S in] AA9069,9531 



(4) 



(5) 



where t = T o /10 4 K. We used the value t = 1.0 for all 
H 11 regions rather than acquiring T from various sources, 
since any variation of t in the range 0.7-1.4 will only re- 
sult in a change in log 77 of 0.1 dex (VP), less than the 
observational errors, and the few regions that have T e de- 
termined in the literature (Webster & Smith 1983; Pagel 
et al. 1979) are all in the range 8000-11000 K. 

In Fig. 6 we plot the deduced 77-parameter for the 
H 11 regions in NGC 300 versus fractional isophotal ra- 
dius which can be compared to Fig. 5 of ZKH. Our data 
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0.20 0-30 

P/Po 



Fig. 6. The softness parameter r\ plotted vs. the fractional 
isophotal radius p/po- The dashed line represents the best lin- 
ear fit. 



points fall in the bulk of their measurements for seven 
spiral galaxies between p/p a = 0-0.5, thus falling below 
the original fit to M33 and M101 by VP. Given the un- 
certainties in all works discussed here this might not be 
significant, though. We confirm for NGC 300 the general 
trend of decreasing log?? with increasing radius. 

VP interpret a correlation of log r] with oxygen abun- 
dance in the sense that at low oxygen abundance stars 
tend to be hotter, as a result of an IMF depending on 
abundance, possibly with some scatter related to evo- 
lution. This implies that for galaxies with steep abun- 
dance gradient one expect to find gradients in and 
likewise in 77, while low-abundance galaxies with no des- 
tinctive gradient only will show a scatter due to mixing of 
evolutinary stages. When compared to the plot by ZKH 
NGC 300 indeed seems to have an intermediate slope of 
the log Ty-gradient resembling maybe that of NGC 2903 
and NGC 4559, which have similar oxygen gradients and 
shallower than the fit for M33 and M101, which also have 
the steepest oxygen gradients. 

The ionization parameter may be calculated from the 
ratio of the [S 11] and [S ill] lines through the relationship 
given by Diaz et al. (1991) 



, ™ , ( [S 11] AA6717,6731 \ „ 
logu = —1.68 log ^ ^ jjjj AA9069, 9531 / ~ 



(6) 



valid for log([S 11]/ [S in]) ~ -0.5 - 1.0. The values of log u 
fall in the range —1.6 to — 3.3 with most of the regions 
scattered around —2.5, thus justifying the use of Eq. (6). 
No trend with isophotal radius was seen as was also the 
case for the seven spiral galaxies where u has been inves- 
tigated by ZKH. 



6. Conclusion 

The observations of NGC 300 have demonstrated the ad- 
vantage of the MOS facility for spectroscopy of extragalac- 
tic Hn regions, where 15-20 regions distributed over a 
wide field can be included in only two exposures. In this 
setup the individual slitlcts can be made sufficiently long 
to enable a good sampling of the local sky background, 
which is necessary to obtain precise spectra at all wave- 
lengths in the range 3650-10000 A. 

We have successfully applied a new sulphur abundance 
determination scheme based on the strong [S 11] and [S in] 
emission lines to H II regions in the spiral galaxy NGC 300 
in spite of the very poor resolution attainable at the time 
of the observations. The regions exhibit a clear trend of de- 
cresing sulphur abundance with increasing galactoccntric 
distance with a gradient of —1.25 ± 0.38dex/p o , as well 
as a distinct oxygen abundance gradient. It is assuring 
that the latter is consistent with previous work (Table 4) . 
Quite low S/O-values for NGC 300 are derived in support 
of the decrease found by Diaz et al. (1991) at high metal- 
licity but any certain conclusions on this issue must await 
more data from a larger sample of galaxies. The radiation 
softness parameter seems to have a gradient with a slope 
equal to galaxies of similar oxygen abundance gradients. 

The sulphur abundance gradient found in this project 
could be improved by including regions with p/p D > 0.5, 
preferably at distances at or beyond p . With the large 
field size and sensitivity of the CCD's now used, multi- 
slit spectroscopy will undoubtedly be a powerful tool for 
abundance studies in spiral galaxies. 

In the future we will apply our empirical sulphur abun- 
dance method to Hn regions in other spirals. Hopefully 
further detailed studies of sulphur abundance in bright 
H 11 regions will lead to an improvement of the calibration 
of the suggested empirical sulphur abundance indicator, 
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